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Electromagnetism becomes a nonlinear theory having (effective) photon- 
photon interactions due at least to electron-positron fluctuations in the 
vacuum. We discuss the consequences of the nonlinearity for the force felt 
by a charge probe particle, and compare the impact of Euler-Kockel QED 
effective nonlinearity to the possibility of Born-Infeld-type nonlinearity. 
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1. Lorentz force and Quantum Vacuum as polarizable medium 

Given the weakness of gravity, a natural question is, if nonlinear elec- 
tromagnetic interactions can on stellar scales lead to transport of matter in 
ways differing from the Lorentz force expectation. This question is espe- 
cially pertinent for stellar objects having extreme magnetic fields (magne- 
tars), and we investigate here in a covariant formulation additional forces 
induced in extreme magnetic fields. 

The Lorentz force predicts that a charged particle moving in a magnetic 
field of arbitrary strength will experience a force normal to both the field 
direction and the direction of motion. This pivotal property decides the 
fate of charged particles in the magnetic field-filled Universe. Even a very 
small field suffices to determine the particle dynamics when acting over a 
distance large compared to the natural microscopic particle scale, and does 
so as a function of their velocity. Electric fields play a subdominant role on 
macroscopic scale since their presence requires a separation or imbalance of 
charge which normally cannot be maintained for a long time or/and in a 
large volume. 
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Fig. 1. Feynman diagram expressing effective field- field interaction going tlirough 
electron-positron fluctuation. 

The Lorentz force has not been well tested on a macroscopic scale. It 
could be modified even on planetary scale, considering the material con- 
stants of a planetary space plasma, the dielectric polarizibility e and mag- 
netic permeability The question we wish to address here is if, and under 
what conditions analogue physical effects originate in the empty space vac- 
uum polarizibility due to vacuum fluctuations and/or a natural nonlinearity 
of electromagnetism. These effects are normally studied in atomic environ- 
ments, and in comparison to the dominant Coulomb force are exceedingly 
small. Precision atomic experiments in part compensates this smallness. 
However as we will argue here the astrophysical environment may offer an 
alternate approach to this physics domain of nonlinear electromagnetism. 

The most famous example of modification of the classical electromag- 
netic theory is the anomalous (since absent unless quantum fluctuations are 
considered) light-light scattering phenomenon. Early in the development 
of quantum electrodynamics, Euler and Kockel [T] (and Heisenberg and 
Euler [2], see below) recognized that electron-positron fluctuations in the 
vacuum would generate an effective self-interaction for the electromagnetic 
field. This effect is visualized in the language of Feynman diagrams as the 
four electromagnetic field legs attached to an electron loop, see Figure [TJ 

One may also view this phenomenon as the presence of the external field 
F'^'^ polarizing the electron fluctuations in the vacuum. This polarization 
of vacuum fluctuations is analogous to the modification of electromagnetic 
properties in a material body, e.g. plasma, but with several important 
differences of fundamental origin: 

1. The QED vacuum is Lorentz invariant, whereas a material body has 
a preferred reference frame in which the medium is at rest. This 
means that, among other things, the framework discussed here and 
the physics it represents provide a consistent Lorentz invariant theory 
of nonlinear electromagnetism with sources. 

2. In QED, the linear polarization 

F^^'^F^, ^ (1 + (a/Svr) ln(-gVm2))F'^-F^, 
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is divergent and is absorbed in the renormalization of charge. In a 
material body, the hnear polarization E ^ D = (1+ eo)-^' is typically 
finite. 

3. The (consequences of QED) polarization we discuss arise from the 
nonlinear response, such as visualized in FiglH and can be presented 
in terms of a dielectric tensor e that is necessarily a covariant function 
of the applied quasi-constant field. 

Heisenberg and Euler showed how to calculate the polarization effect to 
all orders in the applied (external) field [,2j. In so doing, they obtained the 
first non-perturbative result in quantum field theory, which has been stud- 
ied and generalized extensively in the intervening years; for a survey of this 
development we refer the reader to the recent review [3] . The outcome is an 
effective potential Ves for the electromagnetic field, providing the paradig- 
matic example of effective field theory, in this case arising by integrating 
out a "heavy" electron in the low energy [uj <^ rrie) limit of large scale (in 
space and time) quasi constant fields. Their result contains the Euler and 
Kockel result seen in Fig{T] at the lowest finite non-vanishing order (after 
renormalization) expansion in the fine structure constant a. 

However, the full Heisenberg-Euler nonlinearity of electromagnetism is 
important, considering certain compact stellar objects, 'magnetars', which 
are now widely expected to be harboring magnetic fields of equal or greater 
magnitude to the so-called "critical" strength of QED [U [5] 

|4| = = 4.41 10^3 G. (1) 

en 

Near this field strength, all higher order terms in the effective potential 
become of the same order as is the lowest Euler-Kockel term, see figure [TJ 

We will also explore in this report the effective force in a natural nonlin- 
ear extension of Maxwell linear electromagnetism, that is the limiting field 
Born-Infeld electromagnetism [6]. Considering microscopic physics, we show 
that new constraints on the Born-Infeld theory require first understanding 
the dominant QED (Heisenberg-Euler) effects at the percent level. 

In the following, we will recall briefly the framework of classical nonlinear 
electrodynamics and how to derive the dynamics of charged particles and 
fields. We will discuss why it is necessary to consider energy-momentum 
conservation explicitly in deriving the effective electromagnetic force and 
obtain a general formula for the leading terms in the force valid in the leading 
order in the expansion of the nonlinear effective potential. This corrects our 
earlier presentation in theoretical detail [7j, though the magnitude of the 
effects remains. 
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2. Classical Dynamics 

The Lorentz equation of motion of particles, and Maxwell field equations, 
are derived by minimizing the total action. 

/tot = j d^xCtot, Ctot = -S + V^s + fAf, + m^iu^u^ - 1). (2) 

Here the Maxwell electromagnetic action —S = {E'^—B'^)/2 is supplemented 
by an effective potential, such as the Heisenberg-Euler Veg. The Maxwell 
electromagnetic field equations follow from the Euler-Lagrange equations, 
varying Eq. ([2]) with respect to A^. 

Note that the Lorentz force emerges as consequence of variation of the 
path a charged particle takes in the field-filled space. However, a more uni- 
fied view of paths with fields arises when we consider the Lorentz force 
as emerging from conservation of energy-momentum tensor. The Feynman 
path integral formulation of QED offers yet greater unification of both dy- 
namical elements (fields and paths) but this discussion goes well beyond the 
scope of this report. 

At first, we will not specify what V^g is, in order to provide a general 
derivation and discussion of the consequences of the presence of any V^s- 
The third term in Eq. ([2]) is the gauge invariant minimal coupling of the elec- 
tromagnetic potential to the current, and the last term encodes the charged 
particles' (inertial) mass m and four- velocity = 7(1, u). Note that at this 
point we do not differentiate between electromagnetic and mechanical mass 
of a particle: there is inertia in the electromagnetic field that accompanies a 
charged particle and a redefinition of mass is required. This will be achieved 
below in Eq. (|10p where we consider the force related to the nonlinearity of 
the field. 

The effective electromagnetic Lagrangian 

C = -S + V,s (3) 
can depend only on the scalar and pseudoscalar Lorentz invariants 

Varying Eq. ([2]) with respect to A^ then yields 

d,H>^''=f (5) 
where H^'^ defined as the displacement field tensor (see for example §6 of [8]) 

dF^i, dS dV 
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When Ves —> 0, note that H^'^ —> F'^", and Eq. ([5]) becomes the standard 
Maxwell equation with source, dfj,F^'^ = j'^ . 

Current conservation means the field tensor is the 4-dimensional curl of 
a vector potential and implies the homogeneous Maxwell equation 

a^F'^^ = 0. (7) 



3. Charged Particles in External Fields 

The dynamics of the charged particle plus field system can be derived 
from the statement of joint electromagnetic and matter energy- momentum 
conservation, 

MTo'L+TZtter)=0 (8) 

The matter (last) part of the action Eq. ([2]) defines Tj^atter with the property 

d>.TZtter = u^^d,{mun = dp^^/dr. (9) 

At this point, one might be tempted to identify —d^Te.m. as the force 
applied to the charged probe particle. However, this calculation yields the 
standard Lorentz force j^F^'^, even in the presence of nonlinear electro- 
magnetism (see Appendix [X]). The correct approach to recognize the force 
is inherent in the work of Born-Infeld [6J , which considers the field energy- 
momentum of individual particles to be identified with their electromagnetic 
inertia, and hence included in the definition of inertial mass m of the parti- 
cle sourcing the field. We extend this also to the case of vacuum fiuctuation 
nonlinearity. 

By subtracting the energy-momentum of the separate probe particle and 
external field components as they appear in isolation, only the nonlinear 
interaction part of the electromagnetic field tensor is retained as sourcing 
the particle-particle force, see section 4 in Ref. [12j . 



TZ = Ttm.-T^'' -Tr- (10) 

The electromagnetic force is exhibited as the divergence of interaction energy- 
momentum TJ^t • This procedure will show the usual Lorentz force, as well 
as further contributions due to the nonlinear interaction of the charged par- 
ticle's electromagnetic field with the external field, see Appendix iBl 

For a general nonlinear electromagnetic theory the energy-momentum 
tensor can be written [9j (see also Eq. (5) of |10j) 

T^m. = eTZ^ + 9'''T/A, = g^'^S - F^^F\ (11) 
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Fig. 2. Diagram expressing contribution to the effective force arising from light- 
by-light scattering in QED. In the Born-Infeld theory, the loop is to be seen as 
coUapsed to a point. 



The deviation from the Maxwell energy-momentum tensor T^^^ is described 
by the two Lorentz-scalar functions 



d{-S + Feff) _ dVes 
dS dS ' 



(12) 



r^T; = -4(v;,-5^-7'^) (13) 

The energy- momentum trace T must arise from T4ff, because the energy- 
momentum tensor of the Maxwell theory T^^^ is traceless, and it must 
start at order e^, being related to the presence of (effective) field-field inter- 
actions [TI], such as displayed in Fig. [T] 

For the case of an external magnetic field, the interaction energy-mo- 
mentum T'l^ is calculated to leading order in Appendix [Bj This depends 
on three powers of the external field and is order e^, which one can see in 
Fig. [1] by counting the vertices, each of which comes with a power of e0 
In the case of QED nonlinearity we are thus finding the effect represented 
schematically in Fig. [2l 

Taking the divergence of Eq. ()B.7p . 



Ik 



-d,TZ = -j.Fr + Sr (14) 



For the second equality, the Lorentz force (density) is obtained from the 
divergence of the Maxwell-Lorentz interaction Eq. (jB.Sp and separated. 



^ We count powers with the QED Heisenberg-Euler Vcs in mind, and we will see that 
the count is simply decreased by 2 for Born-Infeld theory IS! due to the absence of 
a/n arising from the fermion loop. 
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The remainder is the sought-for modification 

Sf"^ %KxHfd,TZ^_, + TZ,A{%nxH;'') (15) 

where the Maxwell energy-momentum tensor T^^^ ^ and partial derivatives 
de/dS,dT/dS are evaluated at the external field. 

4. Born-Infeld — Euler-Kockel Comparison 

The action of Born-Infeld electrodynamics is 



Cbi = {£c? - {£c?y^TT2S/£f^ry£! (le) 

Born selected originally the limiting electric field strength, £c, such that the 
inertia of the electromagnetic field of an electron was equal to the electron 
mass. Appearance of other fundamental charged particles, the quantum- 
field theory path to mass renormalization, and further mechanisms (such 
as Higgs) to generate mechanical particle mass reveal that Sc = M^/e is a 
free parameter, constrained by experiment [13] . Limitations of the possible 
range of values of the associated mass scale M will be discussed below in 
comparison with QED. 

The expression for the 4-force Eq. (|15p is given in terms of e and T, 
which we can obtain analytically from Eq. ()16p . The coefficient functions 
appearing in Eq. ([22]) are 



^ = -(1 + 2S/S! - V'/S^r'/'S-' = -e'£-' (17) 



M4 MS M12 2 

4dS ^ '^'^^ 



e\S/£l-VVSt) (18) 



given also in weak- field expansion in the second line. 

The Euler-Kockel effective potential, which is the first term of Euler- 
Heisenberg power series in /m^ is given by 
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The common single power of a/vr is due to the calculation being a one- loop 
evaluation of electron fluctuations. For the coefficient functions, one finds 

de ^ 4 a 32 a 

as ^~ 45^^^ 105^^ ^ ^ 



---4^ - —-4 (245^ + 137^^) + - (21) 
4 55 45 vr m4 315 vr ^ ' ^ ' 



Let us compare Born Infeld results Eq. (jl7p and Eq. (|18p to Euler-Kockel 
result Eq. (pOj) and Eq. (j2T]) : we see that the QED effects are apparently 
more strongly suppressed because they arise from quantum corrections, sig- 
naled by the presence of a/vr. With the smallness of a/vr, the Born- Infeld 
coefficients are ~ 5000 (mg/M^) times the QED coefficients and clearly this 
is excluded by precision QED tests. In fact the latest published constraint 
on Born-Infeld M > 60 MeV [13], meaning m^/M^ < 6 10^^, which assures 
that QED dominates Born-Infeld corrections. 

5. Physical efTects of the nonlinear force 

The force derived from the nonlinear action of QED or/and Born-Infeld 
type theory will originate in an increase in field energy when superposing 
fields. Thus in general, the force should act in such a way as to screen 
and reduce the total field strength. In the case of a magnetic-only external 
field, this means the sign of the force depends on Be ■ Hp displaying the 
relative orientation of the external field Bg and probe magnetic field Hp. 
For example, when the fields are aligned, B^ ■ Hp > and the force acts to 
expel the probe from this domain. On the other hand, when the fields are 
anti-aligned, the force acts to pull probe deeper into the strong field domain 
to increase the screening. 

The effect is manifested in the 4th (time-like) component of the force 
determines the change in energy of the probe current and is related the 
3- force given shortly by the covariant requirement u^dp^^/dr = 0. In the 
0-component of Eq. (jl5|) we take the external electric field to be negligible, 
so the first term in Eq. (fT5|) does not appear. With F^^Hp'^ = 2Bc ■ Hp, we 
obtain 

6f = 5e dt (^2^e • Hp^ - dt (^1^24 • 4) (22) 

de -i -± \ „de -t 



ASedt { -^B, ■Hp\+ 2— Be ■ HpdtSe. (23) 



To obtain the second line, we use the identity 



1 dT ^de _de , _ . , 
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which follows from the respective definitions Eqs. ()12p & (|13p . Again, since 
the external electric field is vanishing, the second term in this identity pro- 
portional to V = E ■ B does not appear. 
The 3-force corresponding to Eq. (|22|) is 

- - ae - - - 
6f = 2B, ■ Hp—B, X (V X B,) (25) 

+ (|5ePV - 4(4 • V)) (^2Be ■ i?p^) + V (^2Be ' ^P^f^) 
= 2Be • ^p^^e X (V X B,) - 2B, ■ Hp — VSe (26) 




The second equality is obtained again using Eq. (j24p . The first term in 
Eq. (|26p corresponds to the first term in Eq. (jlSp and as seen in consideration 
of (J/'^ does no work on the probe current. It remains to be shown that the 
remaining two terms in Eq. (|26p are physically relevant. 

We consider two competing forces in an astrophysical plasma near a 
compact star: the gravity of the star and the internal energy (pressure) 
of the plasma. For order of magnitude estimates of the force Eq. ()26p . we 
can use the Newtonian approximation for the star's gravity and replace the 
gradients by l/L where L is the compact star length scale, i.e. its radius 
~ 10 km. Then, the ratio of radial forces is 



(5/1 de\Bo\^\Hp\ L2 _ ia {B^l'^ \Bc\\Hp\ L 



r-/grav| 9S L GMp 457r|4|2 p GM 



(27) 



where p is the mass density of the plasma p ~ rriplO^^ cm~^ for a neutral 
astrophysical plasma. G is the Newton gravitational constant and M the 
mass of the star, so that L/GM ~ 5 — 10 for a compact star. We have 
replaced de/dS by its QED value so that \Bc\ is the QED "critical" field 
strength Eq. ([1]). \Hp\ > a few mG according to estimates and models of 
accreting plasmas [15j . Putting in the numbers shows that the nonlinear 
electromagnetic force dominates gravity for |Se|/|-Bc| ^ 10~^, confirming 
the previous analysis [7]. We estimate the pressure P in the plasma from 
the ideal gas law and compare it to the interaction energy density, which is 
essentially the force derived here without the gradients. Thus, the ratio is 

- dS nT ^^^> 
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where n is the number density 10^^ cm~'^) and T is the temperature, 
suggested by observations to be of order 10 — 100 eV [16]. Since this differs 
from the previous estimate Eq. (j27p only by the factors rrip/T ~ 10^ and 
GM/L, the immediate conclusion is the nonlinearity of the electromagnetic 
interaction may be more important than the plasma dynamics near to the 
surface of the star. Note of course that a dipole field falls off rapidly as 
I Be I ~ and these nonlinear field effects will be subdominant farther 
from the surface, say at 200 times the stellar radius where |-BeP has fallen 
to 10~^ its value at the surface. 

6. Conclusions 

After reviewing the framework of nonlinear electromagnetism, we showed 
here how to calculate the force arising from the nonlinear field-field inter- 
action, shown schematically in Fig. [21 The procedure involves separately 
identifying the polarization of energy-momentum tensor related to presence 
of two components, here external field and particle field, so that we can 
account for the electromagnetic mass of particles. 

With only magnetic fields and Lorentz-Maxwell electrodynamics, the en- 
ergy of charged particles near stellar objects is determined solely by gravity. 
We have shown that once Euler and Kockel [T] (effective) QED nonlinear 
electromagnetism is accounted for, magnetic fields can do work on charged 
particles and this effect beats out gravity for fields that are quite strong 
yet still very far from critical. This insight may have considerable physics 
impact, considering that there is no a priori limit to the magnetic field that 
a ferromagnetic star can source. It is for this reason that exploration of the 
physical consequences arising in ultra strong magnetic fields remains a topic 
of current intense discussion relevant to extreme astrophysical events |17] 
and in various areas of elementary matter physics [TSl 113 [Ml 1211 [221 [23] ■ 

To close we note that the current report must be seen as a first step on 
the way to understand the classical charged particle dynamics in presence 
of ultra strong external magnetic fields and gravity. Beyond the effects 
we considered here we further expect radiation reaction to be relevant in 
proximity of the critical field strength where acceleration is so strong that 
the radiation field impacts the source dynamics. 

Appendix A 

Finding the Lorentz force in NLEM 

In this appendix, we calculate the divergence of the electromagnetic energy- 
momentum tensor Eq. (|lip . 

d.T^:^, = {d,e)TZ^ + ed.TZ^ + d,g^'^T/A (A.l) 
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First, we write out the first two terms 

dC 



{d,e) = [-d,^J {g^'S - F^^F^,) (A.2) 
d,TZ^ = d,{g^^S-F^^^F^,) (A.3) 

with C = —S + 14ff a shorthand for the total electromagnetic Lagrangian. 
Then for the divergence of the trace 7~, we use the form given in Eq. (jl3p . 



W4 = + (a,§)5+(5,|^)p (A.4) 

simplified by virtue of dC = {dC/dS)dS + {dC/dV)dV. Note that it does 
not matter whether one writes C or 14ff in the variation because any terms 
linear in iS manifestly cancel between the first and second terms of Eq. . 
The first term in Eq. (|A.4p then cancels with the first term in Eq. ()A.2p . 
Next, we observe that 

- gf^'V = Ff'^F^" (A.5) 

and find for Eq. ()A.ip 



_ 1 pA*" _ _^ Cn^^5 - F'"' F"" ) + af"^ i B — 
'BS " BS^^^ ^ J + 9 \ci^Q^ 



(A.6) 



'BS 



+ ^{F^''d,F\-g^'^B,S) 



Now using the definition of the displacement tensor Eq. ([6|), we recognize 
the first three terms in brackets as an expanded expression of the divergence 
of —H^'^. Note that Bf^F^'^ vanishes according to the homogeneous Maxwell 
equation Eq. ([7|). From Eq. ([7]) one can also show that 

B^S = ^F^%F,x = -^F'^^B^Fx^ + BxF^,) = -F^^B^Fx^, (A.8) 

which means the final two terms in Eq. (|A.7p cancel. Thus, we obtain 

5^7ri=j^F^^ (A.9) 

A slightly different proof, starting from a different but equivalent form of 
Te'^m. can be found in §6 of [S]. 
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Appendix B 

Lorentz force Correction derivation 

Considering the external field to provide the dominant part of the electro- 
magnetic energy-momentum, we expand in the displacement tensor of the 
probe particle 




L'-f e.m. 



(1) dH^p 



Hf (B.2a) 



(2) dH'^^dH^^ 



-Hp 



T^" (B.2b) 



with the subscript e reminding that derivatives are evaluated for the external 
field. The zeroth order term is just T(f , which is subtracted, and the 
energy-momentum of the probe particle is found at second order in T^g)" 
subtracted. The expansion is in terms of the probe particle displacement 
tensor, because it is defined by the Maxwell equation with source Eq. ([5]). 

The derivatives with respect to H^'^ are related to derivatives with re- 
spect to Ff"" by 

obtained from inverting Eq. ([6]). The 1/2 is a symmetry factor arising 
from the equivalence under permutations of antisymmetric indices. The 
two terms shown correspond to expanding up to order as needed in the 
present approximation. Since T^^ is quadratic in F^'^, expanding Eq. ()B.3p 
in powers of e is the primary expansion involved in obtaining the result here. 

Using Eq. ([IT]) , the calculation of the tensor derivatives in Eq. ()B.2p is 
split into pieces: 



3^f«A Qppa QH^^X Qppa 

de _ de OFP" dT _ dT dFP" 

QJJkX - Qppa QJJkX QJJkX " Qppa QH'^X ^ ■> 

Because V^s and so e^T depend only on S,V, derivatives with respect to 
F^^ can be split into partial derivatives 

j^F,, + —F,, (B.6) 



dFP'' dS ^ dV 
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We see then that de / dS ,dT / dS ,de/ dV ^dT / dV are ah order e^, and the 
next-to-leading-order term in dF/dH Eq. ()B.3P can be dropped. 

Looking to the study of astrophysical magnetic fields, we specialize to 
the case E ■ B = V = Q for the external field. Because V^g must be par- 
ity invariant, it contains only even powers of and therefore the partial 
derivatives de/dV and dT /dV vanish being proportional to V. In this case 
then 

Bf 1 Fit 

T-Z = - TZ.,eQ^K,Hf + 9'''lQ^l'a,Hf (B.7) 

Here ^ is the Maxwell energy-momentum tensor, compare Eq. ()lip . for 
the external field and the (linear) Maxwell-Lorentz interaction is separated: 

Ttp' = -{FtH;^ + K^H^'')g.x + g'^'lK^Hf (b.s) 

By going to the rest frame of the probe particle, it easy to show that this 
term produces the Lorentz force 

d^T^;=j^Fr (B.9) 

Taking the negative divergence of Eq. (|B.7p yields Eq. ()15p . 

Now, studying the second order term T^g) it turns out there are terms 

of order e^. Most can be seen to have form of the electromagnetic energy- 
momentum of the probe Tp'^ (in its nonlinear form Eq. (jlip ). One term 
has the form of a scalar modification to the Lorentz force, which is seen by 
differentiating e and each once. However, we will consider the probe 
field to be much weaker than the external field, in which case these terms are 
all smaller than those calculated by a factor |//p|/|i?e| <C 1. Incorporating 
these terms would provide the complete result for the 4- force at order e^. 
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